On the Canonical Formalism for a Higher-Curvature Gravity 



Yasuo Ezawa 1 , Masahiro Kajihara 1 , Masahiko Kiminami 1 , Jiro Soda 2 and Tadashi Yano 3 

1 Department of Physics, Ehime University 
Mtsuyama 790 Japan 
2 Department of Fundamental Sciences, FIHS Kyoto University, Kyoto 606 Japan 
3 Department of Electrical and Electronic Engineering, Ehime University 

Matsuyama, 790 Japan 



Following to the method of Buchbinder and Lyahovich, we carry out a canonical for- 
malism for a higher-curvature gravity theory in which the Lagrangian density £ is given in 



terms of a function of the scalar curvature R as C = J—detg fiV f(R). The local Hamilto- 

nian is obtained by a canonical transformation which interchanges a pair of the generalized 
coordinate and their canonical momentum coming from the time derivative of the metric. 



§1. Introduction 

The higher curvature gravity (HCG) theory is derived naturally from the superstring 
theory as its low energy effective theory. HCG is one of the modifications of the Einstein 
gravity. For macroscopic phenomena, the Einstein gravity works sufficiently well so that 
the modification is negligible. The modification is, however, expected to play an important 
role in microscopic, i.e. Planck scale phenomena. In most of the unified theories of fun- 
damental physics, such as superstring , supergravity theory and Kaluza-Klein theories, etc. 
extra-dimensional spaces are necessary. The extra-dimensional spaces are expected to have 
contracted to Planck scale during inflation. Then there arises the problem of the stabilization 
of the extra-dimensional spaces which has been investigated by many authors [|T], 0, |^, |4|, [5| . 
In the Einstein gravity, it seems to be difficult to stabilize the extra-dimensional spaces after 
inflation |], |7|]. So it is natural to ask whether these spaces are stabilized in HCG. 

If HCG plays a crucial role in the stabilization, the mechanism would be quantum me- 
chanical. Then we are lead to the canonical formalism of HCG which has been investigated 
by Boulware|| and by Buchbinder & Lyahovich (BL)||. The latter authors generalize the 
formalism of the higher-derivative theory by Ostrogradski to an appropriate form for HCG. 
They modify the Lagrangian density, (2.1) given below, by using the Lagrange multiplier 
method which takes into account the definitions of new coordinates and their relations to 
the time derivatives of the original coordinates. The advantages of the method are that, 
once the modified Lagrangian density is constructed, all the equations are derived by the 
variational principle and the consistency conditions, and that it is not necessary that the 
orders of the highest time derivatives of the coordinates are the same. Thus the procedure 
is straightforward although tedious. They found that, if the additional term in the modified 
Lagrangian contains only the square of the scalar curvature, there exist second class con- 
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straints. However, they have not obtained the Lagrange multipliers and the final form of the 
Hamiltonian constraint. 

In the formalism of Boulware, the Lagrangian density is restricted to the case in which 
it contains terms constructed from the Ricci tensor. The generalization to the Lagrangian 
density given by (2.1) is not straightforward [IT, 12]. 



In this work we will complete the formalism of BL. We reduce the phase space by using 
the second class constraints. We have only one remaining coordinate arising from the time 
derivative of the metric||. The last step is the canonical transformation which interchanges 
the new generalized coordinate and its conjugate mementum. This transformation avoids the 
non-locality when we proceed to the quantum theory. The Hamiltonian thus obtained has 
more natural dependence on the coordinate and momentum than before the transformation. 
We will also comment on the relation to the so called equivalence theorem [|I3], F4 well known 
in classical HCG . 

In section 2, we develope the canonical formalism for the HCG. In section 3, the formal- 
ism is applied to the Robertson- Walker spacetime. Section 4 is devoted to the summary and 
discussion. The problem of stabilization of the extra-dimensional spaces will be investigated 
in a separate paper. 

§2. Canonical Formalism 

In this section we present a canonical formalism of a HCG whose Lagrangian density is 
given by 



C = ^-detg,J(R) (2.1) 

where f(R) is a function of the scalar curvature R. In terms of the ADM variables ]10[, the 
three dimensional metric hij = ^ , the lapse function iV and the shift vector N l , the 
scalar curvature R is expressed in the following d + 1 form: 

R = 2g mn £ n K mn + K 2 - 3K mn K mn + R - 2N~ 1 AN. (2.2) 

Here K k i is the extrinsic curvature 

1 1 

K k i = -£ n h k i = 2jy(d h k i - iVfc-j - N l;k ), (2.3) 
n M is the normal to the surface t = constant 

n' t = N- 1 (l,-N m ). (2.4) 
and C n is the Lie derivative with respect to the normal n^ 1 

C n K kl = N-\d K kl - N j K kl . }j - N{ k K fl - N{ t K k3 ). (2.5) 



In the above (and hereafter), a bold face letter denotes the trace of a 2nd rank tensor, e.g. 
h k iR kl = R . 

Now following BL, we introduce new generalized coordinates 

Qij = (2.6) 
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and "velocity" variables 



d Q 



d N\ and v = d N. 



(2.7) 



Then there are relations between the new coordinates and the time derivatives of the original 
coordinates 

<U>,i = 2NQ ij + Ni-j + Nj-i. (2.8) 
In terms of the new variables, Lagrangian density (denoted as Cq) is expressed as 



where 
and 



C Q = y /deth ij Nf(R Q ). 
R Q = 2N~ 1 h kl v kl -U Q - 2N-'AN - 2N-\N k Q. k + 2N k ' l Q kl ) 



Uq — 3Q mn Q mn — Q 2 



R. 



(2.9) 

(2.10) 
(2.11) 



Using the method of Lagrange multiplier, modified Lagrangian density £*, which takes into 
account the relations (2.8) and the definitions (2.6) and (2.7), is given as follows: 



C* = C Q +p kl {d h kl - (2NQ kl + N k;l + N l;k )} 

+P kl (d Qki - v H ) + P k (d N k - v k ) + P(d N - v). 
The Hamiltonian density 7i* corresponding to C* is given as usual 
H* = P kl d Q kl + P k d N k + Pd N + p kl d h kl - C* 

= p kl {2NQ kl + N k . tl + N l]k ) + P kl v kl + P k v k + Pv-C c 
The action S to be used in the variational principle is given by 

S = J d d+1 xC* 



(2.12) 



(2.13) 



(2.14) 



The variations of the velocitry variables give the "momenta" P kl , P k and P which are 
expressed asR 

pki _ — 2Vhf(R Q )h kl 



p, 



p 



dC Q _ 

dv 



[2.15) 







where /' = df/dR and h = det hij . The velocity variables in (2.13) should be replaced by 
these momenta. Equation (2.15), however, cannot be solved for all velocity variables, i.e. 

1 In the Ostrogradski's formalism, these equations define the momenta canonically conjugate to the 
highest time derivatives of the generalized coordinates. 
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there are primary constraints. Decomposing the second rank tensor (density) into a traceless 
part (denoted by a tilde) and a trace part 



we obtain the following expressions for these constraints 



(2.16) 



C 



ki 



p kl « o 



C k = P k ^0 
C = P w 



(2.17) 



where /c, / run from 1 to d and ps means that the quantity vanishes due to the constraints. 
The only velocity variable that can be solved from (2.15) is the trace of v k i 



N , 



(P/2dVh) + U Q ] + AN + N k d k Q + 2N k ' l Q 



hi- 



where P , the trace part of P kl , is given as 

P = 2dVhf(R Q ) or FL Q = f'-\P/2dVh) = ^j{P/2dVh) 
Using (2.19), we can express H* as 
H* = C kl v kl + C k v k + Cv 



(2.18) 



(2.19) 



+N 



2pHQkl + ^f^l 2d ^) + lf U Q ~ Vhf(4>(P/2dVh)) } + ^PAiV (2.20) 



+2p kl N k]l + -P(N k d k Q + 2N k ' l Q kl ). 

The canonical Hamiltonian density 7{q is obtained by putting 

H* = n + C kl v kl + C k v k + Cv. 



(2.21) 



The velocity variables on the right hand side cannot be expressed in terms of the momentum 
variables and play the role of the Lagrange multipliers. Hq is written in the form 



Here 



H 



Ho = NTi + N k Ti k + divergent terms. 



2QtiP kl + ^PmP/2dVh) + U Q ] + ^AP - Vhf(^j(P/2dVh) 
1 



(2.22) 



(2.23) 



H k = -2h km p m l + - [Pd k Q - 2{PQ kl f 



The Poisson bracket (PB) form of the canonical equations of motion are derived from 7i* if 
we treat the pairs of the generalized coordinates and the Lagrange multipliers, (P kl ,Q k i), 
(p kl ,h k i), (P,N), (P k ,N k ) , as canonically conjugate pairs. Then it is easily seen that the 
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Poisson brackets(PB's) among the primary constraints vanish. Furthermore PB's of N k 7i k 
and the primary constraints also vanish. This means that the time variation of the primary 
constraints are solely given by N7i. Thus the conditions that the primary constraints should 
be preserved are that the PB's of NH. and the primary constraints should vanish and are 
expressed by the following set of secondary constraints as 



r c = n « o 



(2.24) 



c h _ p ki _ 1 ph ki + LpfQkl _ -h M Q) « 

a 2d v a ' 



The first and the second constraints are interpreted as the "Hamiltonian" constraint and the 
"momentum" constraint, respectively. 

The PB's of the Hamiltonian constraint and the primary constraints also vanishes ow- 
ing to the secondary constraints. The PB's of the momentum constraint and the primary 
constraints vanish, whch is desirable if the former one is the generator of the spacial coor- 
dinate transformations, which is to be seen. The third one gives the momenta canonically 
conjugate to h k \ which is not given by the equation like (2.13) but was introduced as a La- 
grange multiplier. Only non-vanishing PB's among the primary constraints and C kl is the 
following: 

{C kl hc),C ij (x')} PB = — 7 P(h kl h lj - -h kl h ij )5(x - x'). (2.25) 

2d d 

The constraints C, C k are first class and C hl , C kl are second class constraints 0. Thus 
the velocity variables v k i are determined from the condition that the secondary constraints 
C kl ~ should be preserved in time: 



v kl = 2N[ (P x p - < L-lq)Q kl + Q km Qr - - d h kl Q mn Q mn - l -R kl 

-P- l (V k Vi-^h kl A)(NP) 
+Q km N m l + Q lm N™ - (Q k m N l]m + Q™N k . m ) + P~\P N%Q kl 



(2.26) 



+2P- X (p + IpQXiV^ + N kk - ^h kl N j ;j ). 



Apart from divergent terms, we have 



d 



C kl v kl = 2N[ (P x p - ^Jlq)F"Q h + ±P(P-Ip kl ) 



, T}klr\ f) m 1 pkl p 
;kl + " (4kmS4l ~ ~r -KfcZ 



-\-N 



^{{P^kQ m) \n — (P mn Qmk);n} + P(P 1 P mn Qmn);k 



+4{p- 1 ( P +^PQ)P i J ;i 



(2.27) 
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Then we obtain the following form for H* 

H* = C k v k + Cv + NH + N k H k + divergent terms (2.28) 

where 

H = 2p kl Q kl + - Vhfty) + P x (p + ^PQ)P kl Q kl + 2P kl Q k m Qi m 



+ lf(ZQmnQ mn + ^Q 2 ) + P(P^P k %i + ^AP - (P kl R kl + ^PR) 

H k = -2h kmP m l + - d [Pd k Q - 2(PQ kl y l ] + P(p- 1 P m "Q m „) ;fc 



^2.29) 



1 



+2[(F* k Q l m U - (P ml Qmkh\ + 4[P" 1 (p + -?Q)P k 



2d 



It appears from eq.(2.29) that it would be more convenient to choose Qij and Q as generalized 
coordinates than Qij. Since we have 

i^Q;, = i" J Q :j + ^PQ + \(PQ tJ + P ij Q)ki 

corresponding canonical momenta are P % \ P/d = P and p' 1 ^ which is given by 

p 'ij = a + p Q ij + IpijQ = a + p/gij + -hVQ) + IjWQ. (2.30) 

It is easily seen that these variables obey the usual Poisson brackets among the canonical 
pairs of variables. Then H and H k in (2.29) are expressed as 

H = 2p' ij {Qij + ^hijQ) + \H - ^hfW 



ipg^-^LlpQ^AP-I] 



+ lp-ip'piig. . _ <±±± Q pi3Q.. + 2 P i *Q i l Q jl + !>(!> P"), ; - P fe <P M ( 2 - 31 ) 
=-2W m : + PQ ;fc 

+p(p- 1 p" m g mn ) ;fc + 2(p fc ; 4 m - p' m Q fci ) ;m + ^(p-'p'p^);; 

To proceed further, one of the two methods is usually adopted. One is to use only 
independent variables by solving the constraint equations . The other is to use the Dirac 
brackets without solving the constraints, which is sometomes more convenient when we 
cosider symmetry properties. Here we adopt the former method since it is simpler. As to 
the latter one, we only present the list of the Dirac brackets in the appendix. 

Now we reduce the phase space by using the second class constraints, P u ~ and 
Q l i ps 2P~ 1 p' l i . Remaining independent canonical pairs of variables are {hij,p' 1 ^) and 
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(Q,P), in terms of which H and Hk in (2.31) are simplified to 



H = 2P- 1 ( P '% - ip' 2 ) + ^p'Q + ifV - y/hfty) 



d' 



d 



H k 



±ipQ 2 + AP-i P R 

-2h kmP '™ + PQ k 



(2.32) 



In classical theory, (2.32) gives the final form for the Hamiltonian constraint and momen- 
tum constraints. However, when passing to the quantum theory, we meet a difficulty with 
(2.32) in that it contains terms with negative power of P. Such terms lead to non-locality 
in the canonical quantum field theory. Thus we make a canonical transformation which 
interchanges the momenta with the coordinates 



(Q, P) — (Q', P') = (P, -Q) 



(2.33) 



In terms of the new coordinates and momenta, tp is a function of Q' j2\fh and H, H k are 
expressed as 



H = 2Q- 1 (p'% - ip' 2 ) - ^ P 'P + ^ - Vhfty) 

-^QP 2 + AQ-iQR 

H k = -2h kmP >™ - QP k 
where we have dropped the primes denoting the new variables for simplicity. 



(2.34) 



§3. Robertson- Walker spacetime 



In this section we apply the results of the previous section to the Robertson- Walker 
spacetime whose metric is written as 



ds 2 = -N 2 (t)dt 2 + a 2 (t) 
The scalar curvature R is given by 

R = QN~ 2 



dr 2 



1 — kr 2 



+ r 2 (d6 2 + sin 2 6d<j) 2 ) 



(3.1) 



d /d\2 kN 2 Nd 
- + (-) 

a v a y cr TV a 



(3.2) 



The 3-dimensional scalar curvature is R = 6k /a 2 . The extrinsic curvature is given as 



Kij ~ Na Hl] 



(3.3) 



where hij is the 3-dimensional metric tensor and has a form a 2 (t)x (fixed functions of 
r,9,(p) as can be seen from (3.1). 
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The new generalized coordinates Qij can be taken to be K,- L j as in the previous section: 

Qij = Kij = hij. (3.4) 

Obviously there is only one independent component which we take to be the trace part Q : 

Q = 3(a/Na) or = ^Q/iy. (3.5) 

Thus we take Q to be the new generalized coordinate and the corresponding velocity vari- 
ables as 

v = Q. (3.6) 

Then the scalar curvature is expressed as 

R = 2N-\+ A -G? + % = R Q . (3.7) 
3 cr 

Using this Rq , the Lagrangian density is written as 

C Q = ^f(R Q ) (3.8) 

The modified Lagrangian density, which takes into account the relation (3.5) and the defi- 
nitions such as (3.6), is given by 

£* = C Q + P(Q - v) + P{N -v)+ p(a - ^NaQ). (3.9) 

Here P was denoted as P in the previous section and p is related to p' 1 ^ as p = 2hijp hJ J a. 
It is noted that all the fields are homogeneous. The Hamiltonian density corresponding to 
this C* is taken to be 

H* = PQ + PN + pa -£* 



(3.10) 



(3.11) 



= Pv + Pv + \ NapQ - C Q 
Variations with respect to the velocity variables lead to the following equations 

P = ^ = 2Vhf(R Q ) 

P= 6 -^ = 0. 

ov 

From the first equation of (3.11), we have 

v = ]-N[i)(P/2Vh) - ^Q 2 - 6k/a 2 ] (3.12) 

where 

V(P/2v^) = f'-\P/2Vh) = R Q . (3.13) 
The second one of (3.11) gives the primary constraint 

C = P^0 (3.14) 
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Then H* is written as 




(3.15) 



where 



C Q = NVhf(ij(P/6Vh). 
Using (3.12) we obtain the following expression for 7i* 



(3.16) 



H* 



Cv + NH 



(3.17) 



where 



n = \^ ~ t Q2 ~ 6k/a2) + l aQp ~ 



(3.18) 



We can treat each pair of variables 



(a(t),p(t)), (N(t),P(t)), (Q(t),P(t)), 



to be canonically conjugate as in the previous section and calculate the Poisson brackets. 
Preservation of the primary constraints yields a secondary constraint 



This is the so called Hamiltonian constraint. Poisson bracket between C and H* or H 
vanishes, so that we have no additional constraint. 

By the canonical transformation which interchanges the coordinate Q and the momentum 
P, the final form of H is obtained 



where ip = ip(Q/2\/h) . Eq.(3.20) is obtained from (2.34) by using the relations R = Qk/a 2 , 
p' lj = aph' lj I Id and putting d = 3 

§4. Summary and discussion 

Following the method of BL, we presented the canonical formalism of HCG in which 
the Lagrangian density is given by C = y/—gf(R) ■ Using the second class constraints, 
we reduced the phase space which is introduced naturally in the BL formalism. At the 
last step of the procedure, we made a canonical transformation which interchanges the 
coordinate corresponding to the higher derivatives of the original coordinate and its conjugate 
momentum. This transformation is necessary to avoid the non-locality when we proceed to 
the quantum theory. In addition, comparison of the expressions for H , (2.32) with (2.34), 
shows that the dependence of the Hamiltonian on the coordinates and the momenta becomes 
more natural by this transformation. 

The dynamical degrees of freedom of the HCG considered here is just that of the Einstein 
gravity coupled to one scalar field. The expression of H in (2.34), however, does not appear to 
describe such a system. We introduce a scalar field $ and make a conformal transformation 



ow {c,n*}p. B . = -n. 



(3.19) 



n = \W - U> 2 - 6k/a 2 ) - l -apP - y/hfty) 



(3.20) 



$ = k- 1 In (« 2 Q/2v^), hij = (K 2 Q/2Vh) 2/d h i:i . 
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If this is a part of a canonical trasnformation, corresponding canonical momenta II and j? J 
are related to P and p' lj as 



K 

Vh 



1 1 ^ 

5P + & n 



jiftj = exp (2/c$/d) - f-p + — II W J 



Then takes the following form 



H = 



K 

Vh 



P %3 Pij - 



1 



d-1 



P 



k~ 2 R + 



d-1 1 



n 



2k i2 

p 



4rf Vh 1 d(d-l) 



+K- 2 \Jh V - - exp (-k^)F(^) - {exp (2/c$/d) - 1}R 



+^4^V^ exp QKf&IQWd&djQ 



d 

where F = f /2k 2 . This expression is very similar to the one for the Einstein gravity 
coupled with a scalar field. The first two terms are just the Hamiltonian of the Einstein 
gravity. However, the remaining terms are not the one of the scalar field used in the classical 
equivalence argument. So, at this stage, the equivalence is not clear in the quantum theory, 
although there might exist a choice of variables which would show the equivalence explicitly. 

For the Robert son- Walker spacetime, the second class constraint does not appear since 
the phase space is reduced from the beginning. This will make the situation simple when 
we consider the stability of the internal spaces. In this case the second class constraints are 
restricted to those with respect to the internal variables. Consequently the Dirac brackets 
are required only for these variables if necessary. Investigation of this problem would require 
quantum mechanical treatment and will be given in a separate paper. 



Appendix. Dirac brackets 

In discussing the symmerty properties of constrained systems, Dirac bracket formalism 
is sometimes more convenient than the method using only the independent variables. So we 
here list the Dirac brackets among the canonical variables. 

{^(x),p u (x')} OB = (3#<SJ - - d h t3 h kl )5{* - x') 
{/^( x )>M x ')}db = 

{^■(x), P fe V)W = o 

{Q^'(x), P kl (x')}DB = \h ij h kl 5{x - x') 
{Q^'(x),Q w (x')W = 
{^■(x),P«(x')} D b = 
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{h tJ (x),Q kl (x>)} DB = 

{h ij (x.),P k i(-x.')} DB = 2dQ~ 1 (h ik h j i - ^j/i w )<J(x - x') 

{^'(x),Q fe '(x')}DB = -M{h ik y l - V'/i w )<f(x- x') 

a a 



{^(x),P H (x')W 



- P^7i fe - Q^(2p - PQMSj - ^h ij h kl )]6(x - x'). 



The Dirac brackets of hij,p kl and Q,P are given as 



{Q,P} DB = d5(x-x') 

Q}db = 
{^,P}ds = 

{p^QW = -^Q5(x-x') 

a 



{ P v, P} DB = (2P* - -h ij P)5(x - x') 
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